Itinerant and local-moment magnetism in strongly correlated electron systems 
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Detailed analysis of the magnetic properties of the Hubbard model within dynamical mean-field 
theory (DMFT) is presented. Using a RPA-like decoupling of two-particle propagators we derive an 
universal form for susceptibilities, which captures essential aspects of localized and itinerant pictures. 
This expression is shown to be quantitatively valid whenever long-range coherence of particle-hole 
excitations can be neglected, as is the case in large parts of the phase diagram where antiferromag- 
netism is dominant. The applicability of an interpretation in terms of the two archetypical pictures 
of magnetism is investigated for the Hubbard model on a body-centered cubic lattice with additional 
next-nearest neighbor hopping t' . For large values of the Coulomb interaction, local-moment mag- 
netism is found to be dominant, while for weakly interacting band electrons itinerant quasiparticle 
magnetism prevails. In the intermediate regime and for finite t' an re-entrant behavior is discovered, 
where antiferromagnetism only exists in a finite temperature interval. 

PACS numbers: 71.10.Fd,71.45.Gm,75.10.-b,75.30.Kz 
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I. INTRODUCTION 

The magnetic properties of solids are typically de- 
scribed in terms of two archetypical and opposing view- 
points. On the one hand, the picture of weakly interact- 
ing itinerant electron magnetism is usually implemented 
for metallic systems. On the other hand, for small overlap 
between sites or for strong local interactions the valence 
electrons can form localized moments and behave effec- 
tively like spin degrees of freedom of a Heisenberg model 
(for introductory texts, see, for example, Refs. 1 and 2). 

In strongly correlated electron systems, such as cuprate 
high-temperature superconductors 3 or heavy fcrmion 
systems, 4 such a clear distinction is often obscured since 
aspects of both pictures appear. Due to the strong 
Coulomb interaction the electrons are rather localized 
furnishing large magnetic moments. However, at low 
temperatures usually a bandstructure of heavy but itin- 
erant quasiparticles around the Fermi level forms, giv- 
ing rise to the screening of local moments. Central to 
the magnetic properties in these systems is the com- 
petition between quasiparticle-band formation, possibly 
spin-polarized, and ordering tendencies of large local mo- 
ments. 

This competition yields especially interesting and novel 
physics in the case of frustrated systems. 5 In the re- 
cently studied frustrated heavy fermion systems, such as 
LiV204, frustrated local moments dominate at elevated 
temperatures, while at low temperatures strongly corre- 
lated quasiparticles emerge (see, for example, Ref. 6). 
Similarly, a competition between frustrated magnetic 
moments and a low temperature Fermi liquid has been 
proposed as an explanation for the re-entrant Mott- 
transition found in the highly frustrated organic com- 
pound k-(ET) 2 Cu[N(CN) 2 ]C1. 7 In the context of the 
cuprate high-temperature superconductors there is a 
longstanding question whether the magnetic properties 



are best described in the local or itinerant picture, see, 
for example, Ref. 8. Due to the proximity of these sys- 
tems to the Mott-insulating phase, the localized picture 
is often proposed, whereas recent studies also revealed 
the presence of itinerant spin- fluctuations throughout the 
whole Brillouin zone. 9 

In the Hubbard model, all these aspects can be cap- 
tured. Weakly interacting band electrons can be studied 
for small values of the Coulomb repulsion. With increas- 
ing interaction strength the system undergoes a Mott 
metal-insulator transition above which it is best charac- 
terized by an effective Heisenberg spin model. The inter- 
esting region is close to the Mott transition where large 
magnetic moments and itinerant quasiparticles compete. 
Frustration effects can be studied by considering, for ex- 
ample, the influence of a next-nearest neighbor hopping 
in an otherwise bipartite lattice. 

The magnetic properties of the Hubbard model are 
very well studied. Even though the model was introduced 
in order to describe ferromagnctism (FM) this phase is 
restricted to very large interaction strength and depends 
on details of the bandstructure. 1,10 Antiferromagnetism 
(AFM), however, represents a generic phase of the model 
and is found in large regions of phase space. 1 

A very useful method to study the Hubbard model 
in the strongly correlated regime near the Mott transi- 
tion is provided by the dynamical mean field theory 11 
(DMFT). With this method one can access the single- 
particle Green function as well as two-particle quantities 
like susceptibilities. 

We will consider the Bethe-Salpeter equations and de- 
rive coupled equations for dynamic lattice susceptibilities 
which prove equivalent to the usual expression known 
from the DMFT, 12,13 but lend themselves more directly 
to our specific investigation of magnetism. The ultimate 
goal of this approach is to analytical continue these equa- 
tions directly to the real axis and obtain integral equa- 
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tions for functions of real frequency arguments. 

In this work, however, we employ an additional de- 
coupling scheme for the internal Matsubara summations 
which allows us to derive an universal approximation for 
the susceptibility, which unifies localized and itinerant 
aspects of magnetism. We discuss the range of validity 
of this formula and demonstrate its accurateness by re- 
producing the DMFT Neel temperature for the Hubbard 
model on a simple cubic (SC) lattice in three dimensions, 
D = 3. The appeal of this approach lies in its simplic- 
ity compared to the usual Matsubara approach to two- 
particle properties. 14,15 Only functions of real frequencies 
enter and the inversions of Matsubara-space matrices is 
avoided. 

The competition between local-moment and itiner- 
ant quasiparticle magnetism is studied for a three- 
dimensional body-centered cubic (BCC) lattice with an 
additional next-nearest neighbor hopping t'. In situa- 
tions where quasiparticle magnetism is dominant, the 
tendency toward magnetic order should be strongly sup- 
pressed by a finite next-nearest neighbor hopping, since 
the perfect nesting property of the Fermi surface is re- 
moved. On the other hand, local-moment magnetism is 
less sensitive since the geometric frustration induced by 
the next-nearest neighbor hopping reduces the effective 
exchange coupling only quadratically in t' . This behav- 
ior is indeed found for weak and strong Coulomb repul- 
sions, respectively. In the crossover region for interme- 
diate interactions, an interesting re-entrant behavior is 
found where the effective picture to be used depends on 
the temperature of the system. Crucial to the occurrence 
of this behavior is the presence of the frustrating next- 
nearest neighbor hopping t'. The findings of this work 
arc in accord with the common believe that frustration 
is capable of producing novel and interesting phenomena. 

The paper is organized as follows. After a brief intro- 
duction of the model, Sect. II sketches the derivation of 
the Bethe-Salpeter equations for lattice susceptibilities. 
We present two different interpretations of the result- 
ing simple formula for the magnetic susceptibility, one in 
terms of local moments and the other utilizing itinerant 
quasiparticles. In the first part of Sect. Ill we investi- 
gate the validity of our approach and compare the Neel 
temperature for the SC lattice to known results from the 
literature. The second part of this section then focuses 
on the BCC lattice where the magnetic properties are 
investigated in detail. We also include two appendices, 
where more details on the technicalities arc presented. 
The focus is laid on the similarities and differences of the 
the single-particle and two-particle set-up. 



II. MODEL AND METHOD 

In metals the long-ranged Coulomb interaction is usu- 
ally screened and only short-range components have to 
be considered. In order to study the effects of electronic 
correlations we consider the Hubbard model, where only 



the on-site Coulomb repulsion is retained. For the sake of 
simplicity, we assume ionic s-shells without orbital degen- 
eracy and with remaining spin-degeneracy only, although 
the formal developments described below can easily be 
generalized. 

The Hamiltonian is given by 
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e Clavier 



U ^hifflil, (1) 



where the operator c ia {c\ a ) annihilates (creates) an elec- 
tron in a localized Wannicr orbital at lattice site i with 
spin cr, fii a = c ia .c icr is the number operator and e is the 
ionic level position where the chemical potential is al- 
ready taken into account. The electrons can transfer be- 
tween sites i and j with an amplitude Uj , which accounts 
for the itinerancy and band formation of the electrons. 
The last term in Eq. (1) implements the local Coulomb 
repulsion with the interaction matrix element U. 

The structure of the lattice is encoded into the one- 
particle hopping amplitude ty, the Fourier transform of 
which gives the single-particle dispersion relation 
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The noninteracting density of states (DOS) is entirely 
determined by the dispersion relation 
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t k )- 
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Despite its simplicity, the exact solution of the Hub- 
bard Hamiltonian of Eq. (1) is only possible in one spa- 
tial dimension 16 (for a recent book, see, e.g., Ref. 17) and 
in infinite spatial dimensions, D — > oo (with an appro- 
priate rescaling the hopping parameters), 18 by means of 
the DMFT. In this work we will utilize DMFT in order 
to extract an approximate solution for three dimensional 
systems. 

In finite dimensions the major approximation of 
DMFT is to treat all spatially nonlocal correlations in 
a mean-field manner. For the single-particle properties, 
this implies the assumption of a momentum independent 

interaction self-energy, TP (k, z) Dr ^> T YF{z). Then, the 
structure and dimensionality of the lattice enters the lat- 
tice Green function only via the dispersion relation, 



G(k, z) 



1 



-TF{z)-tk 



(4) 



However, the dependence on the (complex) energy 
variable z is not neglected and thus dynamic local corre- 
lations arc fully retained within DMFT. The self-energy 
and the local Green function, 
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can be obtained from an effective single-impurity Ander- 
son model (SIAM) embedded in a self-consistent medium 
characterized by the hybridization function T(z). Once 
the effective SIAM is solved for a given T(z), a new guess 
for the self-energy is obtained by inverting Eq. (6), which 
in turn is used to get G(z) via Eq. (5). Only in this last 
step the lattice structure enters. The self-consistency cy- 
cle is closed by reorganizing Eq. (6) and extracting a new 
guess for T(z). 

The nontrivial part in this cycle is the solution of the 
effective SIAM. But due to its long history, a multi- 
tude of different methods for treating this model exist, 
for example, exact diagonalization, 19 several variations 
of quantum Monte-Carlo schemes, 20 and the numerical 
renormalization group 21 (NRG). In this work, we employ 
the enhanced non-crossing approximation 22,23 (ENCA), 
which has no adjustable parameters and works directly 
on the real frequency axis. 

The tendency toward magnetism is investigated via the 
magnetic susceptibility of the paramagnetic phase, which 
diverges at a second order phase transition. The suscep- 
tibility is given by the nonlocal time-dependent order- 
parameter correlation function 

xr S WHT[M*(r)M;(0)]), (7) 

with the total magnetization operator at lattice site i 



(8) 



In general, the index a denotes orbital and spin quantum 
numbers. With the assumption of s-shells one has espe- 
cially a = <j and j a = — ^HS-cr, where g is the electron 
Lande factor and ^lb the Bohr magneton. 

Equation (7) is already specified to isotropic situa- 
tions where the susceptibility tensor is diagonal and all 



diagonal elements are equal, i.e. x 



mag;a/3 
ij 



Saf: 



= {x, y, z}). For a paramagnetic situation, the 
transverse susceptibility is just given by twice this value, 

mag;± _ 9 v ma S 
Xij Xij ■ 

The wave- vector and frequency dependent susceptibil- 
ity is obtained by the Fourier transform of Eq. (7), 



1 f p 

= y^lalb\(h a )(h b )8iv n .o+ x(q,ivn) \ 

{ L= — J abba J 



(9) 



where v n = =^-n (n G Z) arc bosonic Matsubara frequen- 
cies and fi = the inverse temperature T. We already 
separated the static (iv n = 0) unconnected part propor- 
tional to a product of local occupation numbers (n a ) in 
the hrst term. 

We introduced a matrix notation in orbital space 
for the Fourier transform of the connected two-particle 
Green function, ^(q, iv n ) = {Xa,b;c,d(q,^n)} ( se e Ap- 
pendix B). This is the central quantity of interest, as 



other particle-hole susceptibilities can equally well be cal- 
culated from it. The charge susceptibility, for example, 
is obtained by using different matrix elements j a = — |e| 
in Eq. (8), which amounts to a different sum over the 
matrix elements of the susceptibility matrix in Eq. (9). 
Therefore, we will focus in the following on the suscep- 
tibility matrix xi.1i' LV 'n)i without the specific pre-factors 

7o unless needed, i.e. for their signs. 

In analogy to the self-energy, the particle-hole irre- 
ducible two-particle interaction vertex is assumed to be 
momentum independent within DMFT, 12 ' 13 



A (iwi , k x , iu 2 , k 2 ; iJ 2 ,tf 2 , icj[ , ) 



(10) 



DMFT 



n(iwi, iu> 2 ; iuj 2 , iuj'i). 



As a consequence, sums over crystal momentum only in- 
volve the one-particle propagators leading to simple lo- 
cal expressions for closed loops or a geometric series for 
chains through the lattice. At interaction points crys- 
tal momentum is only conserved on average. A similar 
procedure fails for internal Matsubara sums since all fre- 
quency dependencies are retained within DMFT. There- 
fore, the Bethe-Salpeter equations, as detailed in the ap- 
pendix, have the structure of coupled integral equations 
in Matsubara frequency space. 

The dynamic susceptibility is obtained by summing the 
appropriate two-particle Green function over two internal 
frequencies, 

X(£^«) = 4 E X(£>n|iwi,iw2) e^+^ s . (II) 

The exponential incorporates infinitesimal convergence 
factors S which ensure the correct time ordering of the 
number-operators in the two-particle Green function. 

Frequently, in particular when using QMC as impurity 
solver, the two-particle lattice susceptibility is obtained 
by interpreting Green functions as matrices in Matsubara 
frequency space, x(g, if n |iwi, ia; 2 ) = Xqw, U,^ (indi- 
cated by a hat and the different placement of the external 
variable in our notation). The Bethe-Salpeter equation 
then has the structure of a matrix equation, where the 
internal frequency sums are represented by matrix mul- 
tiplications, 



Xq.; 



-P . 



P.Q,iv n n.iv n Xq,' 



(12) 



The particle-hole propagator P_q, Un is determined by the 

single-particle Green function (cf. Eq. (B29)), and ex- 
plicit two-particle interactions are incorporated via the 
irreducible vertex n, / n , which is a priori unknown and 

hard to calculate directly for strongly correlated systems. 

For a description via effective impurities underlying 
the DMFT-method, an analogous equation can be for- 
mulated, 



y 1oc 



ploc _ p 



n, 



(13) 
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Xiv n is the dynamic local susceptibility of the effective 

impurity model and can be calculated in principle. The 
local particle-hole propagator is given by the momentum 
sum of its lattice counter part, 



ploc _ j_ V P 



(14) 



and can be calculated with knowledge of the local Green 
function (see Eq. (B28)). 

The same irreducible local vertex Tli Vn is assumed for 

the impurity and the lattice model. This makes it possi- 
ble to eliminate it from Eqs. (13) and (12) and the usual 
DMFT result for the lattice susceptibility is obtained 



Xq,- 



-P- L 



Y l0C 



P 



loc 



(15) 



The inversions indicate matrix inversions in the space of 
Matsubara frequencies and orbital/spin space. It is im- 
portant to realize that all quantities entering Eq. (15) can 
be calculated directly from the effective impurity model 
and the lattice Green function. One thus has an explicit 
equation for the lattice susceptibility. 

However, the shortcoming of the approach sketched 
above is twofold. First, the set of Matsubara frequen- 
cies iuj n is infinite as n G Z. The Matsubara matri- 
ces are therefore infinite dimensional, which is of course 
not sustainable in practical calculations and truncations 
have to be introduced. Second, the calculated suscepti- 
bilities are functions of complex Matsubara frequencies 
which have to be numerically continued to the real axis 
which is mathematically an ill-defined problem. There 
exist sophisticated techniques, such as maximum entropy 
methods, 24 but uncertainties and uncontrolled errors re- 
main. 

Another approach utilizes the dynamic density-matrix 
renormalization group technique and a subsequent 
deconvolution. 25 Apart from being very resource consum- 
ing, the deconvolution to obtain data on the real fre- 
quency axis might also introduce artefacts. 

In order to avoid these difficulties we take a different 
route here. We introduce an additional approximation 
and decouple the Matsubara sums in the Bethe-Salpeter 
equations in a manner similar to the random-phase ap- 
proximation (RPA). For technical details see appendix B. 
The major advantage is that the analytic continuation of 
all frequency variables to the real axis, like iv n — > v + iS, 
can be done analytically. As derived in the appendix, the 
result is 26 - 2728 



X(g, v) = 

-£(^)" 1 +x 1 »" 1 +^ ,c W" 



(16) 



The structure of this equation is of course the same as 
Eq. (15), but the important difference is that no Mat- 
subara matrices occur (there are no hats are over the 



quantities) and the external frequency v is a real vari- 
able. The indicated matrix structure and inversions only 
refer to the orbital/spin space. 

In the following we will focus on the magnetic suscep- 
tibility and on simple s-shells. In that case, all quantities 
are matrices in spin-space only. The particle-hole propa- 
gators of the paramagnetic regime are spin-symmetric, 
Pf = Pj. = P, and only two independent compo- 
nents exist for the susceptibilities, Xt,t;T,t = Xl 1-4.4 an< ^ 
Xt44,t = X4,t;T4- The magnetic susceptibility is a lin- 
ear combination of these two components, Xmag(<Z,^) = 
Xt,t;t,t(9) v ) ~ Xt44,t(9' v )i ano - can b e expressed as (ne- 
glecting the prefactor (^ 2 -) 2 ), 



Xmag(9,^) 



1 



1 



1 



P(q,u) P loc (v) x!£ c a>) 



(17) 



All quantities entering this form are scalar functions of a 
real variable. XmagM i s the dynamic local magnetic sus- 
ceptibility of the impurity model, _P loc (z^) and P(q, v) are 
the local and lattice dependent particle-hole propagators, 
which only depend on the fully interacting single-particle 
Green function (see Eqs. (B28) and (B29)). 

This specific form (17) for the susceptibility lends itself 
directly to interpretations in terms of the two archetypi- 
cal physical pictures underlying magnetism: 

(i) The picture of local-moment magnetism is sug- 
gested if Eq. (17) is re- written as 



XmagO?, V) = 



xUlM 



1 - J{q,v 



AmaeV / 



(17a) 



where local spins characterized by Xmag interact via 
a nonlocal and dynamic exchange coupling 



J{q,v) 



1 



1 



P{q,v) P Xoc {v)' 



(18) 



This view is substantiated by approaching the 
atomic limit. Then, the hopping matrix ele- 
ments vanish, tk = 0, and the lattice and the 
local particle-hole propagators become identical, 
P(q, v)\ t=0 = P loc (i/). The susceptibility correctly 
reduces to the fully interacting susceptibility of iso- 
lated ions, 



Xmag (q,f) I 



Xz 



(19) 



For half-filling and U > it is just given by that of 
a free spin 



Xmag (q,v) 



t=o 



5{y) 



(20) 



This result is easily extended to incorporate the 
leading order correction due to the coupling to 
neighboring ions. Expanding up to second order in 
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the hopping around the atomic limit yields the ex- 
change coupling J(q, v) = ^ ^ X)* *k *fc+« + 0(t 3 ). 
The susceptibility then correctly reproduces the the 
mean-field approximation to the Heisenbcrg model, 

Xmag(g^) = T _ l Tc{q) 6(v)+0(t 3 ). (21) 

In case of a simple-cubic lattice with near- 
est hopping, the Tc{q) is given by Tc(q) = 

-TrEiLi C0S fe)- This reproduces the AFM 
Q = 7r(l,l,...) T mean-field transition at the well- 
known critical temperature Tc{Q) = . 

(ii) The interpretation in terms of the itinerant picture 
of magnetism is obtained by re- writing Eq. (17) as 

-P{q,v) 

Xma. s {q, v) = — — -=-— -. (17b) 

1 +T mas {i/)P(q,u) 

Here, propagating electrons characterized by their 
particle-hole susceptibility —P(q, v) interact locally 
with the vertex 

rmagM = " P^Hy)' (22) 

This view is substantiated by the noninteracting 
limit, where the Coulomb interaction matrix el- 
ement vanishes, {7 = 0. Then, Wick's theo- 
rem is applicable and the local susceptibility re- 
duces to the negative local particle-hole propaga- 
tor X m ag(^)| C/= ° = -P loc {v). The susceptibility is 
then indeed that of noninteracting particles on a 
lattice given by the Lindhard-function 

x{g,v)\ U=0 = -P{q,v)\ U= ° (23) 

1 \ ^ f{tk+q) - f{tk) 

N 2^v + iO+ + tk+ g -t k ' 

Improving this by employing a perturbation theory 
in U we arrive at the usual RPA expression. In 
the lowest order the local vertex is nothing but the 
bare Coulomb interaction, r mag (^) = U, and then 

-Plq.u)\ u=0 

The two opposing cases of a weakly interacting elec- 
tron gas and the atomic limit are exactly incorporated 
in the approximate form for the susceptibility, Eq. (17). 
This substantiates the hope to be able to correctly de- 
scribe the regime of intermediate coupling strength and 
especially the transition from itinerant to localized forms 
of magnetism. 

At this point a comment on the advantages and short- 
comings of the final expression, Eq. (17), is in place. The 
major advantage of this formula lies in its simplicity in 



combination with the correct incorporation of the weak 
and strong coupling limit. And as it will be shown in the 
following, the qualitative and quantitative description of 
correlation effects even for intermediate coupling is very 
good. 

In the usual weak-coupling RPA 29 or exten- 
sions thereof, such as the fluctuation exchange 
approximation 30 and the two-particle self-consistent 
approach 31 , the particle- hole propagator is either calcu- 
lated with bare propagators, i.e. without interactions, or 
interaction processes are included only at a very crude 
level. This also applies to auxiliary boson approaches 32 
and equation of motion decoupling schemes. 33 In con- 
trast, the particle-hole propagator employed here is cal- 
culated with the fully interacting single-particle Green 
function obtained from DMFT, which incorporates life- 
time and many-body effects and can already lead to con- 
siderable modifications. Furthermore, in the above men- 
tioned approximation schemes the two-particle interac- 
tion vertices are given by weighted linear combinations 
of bare interaction matrix elements, i.e. by numbers, and 
the frequency dependence of two-particle vertices is usu- 
ally ignored. In the present treatment nontrivial dynamic 
interaction vertices T(v) or J(q,v) are incorporated into 
the susceptibilities. Finite order cumulant expansions in 
the hopping t as done, for example, in Ref. 34 do respect 
the leading frequency dependence of the two-particle ver- 
tex and include leading order nonlocal effects. Lifetime 
broadening and Kondo physics, however, are not cap- 
tured in these treatments. 

Our approach requires a thorough assessment of the 
quality of the RPA-like decoupling of frequency sums, 
which could possibly cause errors. Virtues and limita- 
tions of our result contained in Eq. (17) will be borne 
out by the discussion in the next section. However, 
we can already formulate the expectation, that in situa- 
tions, where coherent two-particle propagations through 
the lattice lead to the build-up of nonlocal correlations 
the decoupling is likely to fail. In the general set-up, the 
coherence is maintained at each interaction vertex due 
to its energy dependence, whereas in the decoupled ver- 
sion averages are taken at each local vertex. Therefore, 
coherence is maintained only between consecutive local 
two-particle vertices, so that long-ranged correlations are 
affected. The most prominent example where such exci- 
tations are crucial is FM, as it is clear from the conditions 
favoring Nagaoka-type FM. 1 ' 35 As will be demonstrated 
below, this leads to a very accurate description of AFM 
phase transitions, while FM transition temperatures are 
overestimated, as in the usual Stoncr theory. 

Another source of inaccuracies stems from the method 
chosen to solve the effective impurity model. In most 
cases, the drawbacks of the impurity solver pose the 
strongest limitations on the accuracy and validity of the 
calculated Green functions. 
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III. RESULTS 

The following investigation of the Hubbard model on 
two different lattices uses the DMFT scheme, applied to 
various one- and two particle properties furnishing, e.g., 
spectral properties, local moments, susceptibilities, effec- 
tive interactions and transition temperatures. We con- 
centrate on the approach from the paramagnetic regime 
in order to identify phase transitions via diverging mag- 
netic susceptibilities. We set <?^b = h = c = = 1 
and use the noninteracting half bandwidth W as unit of 
energy. 

For the results presented in this Section, we used the 
ENCA 22 as impurity solver in order to calculated the lo- 
cal self-energy and susceptibility. The three-dimensional 
momentum integrations of the lattice particlc-hole prop- 
agator were performed using the CUBA-library. 36 



A. Simple-cubic lattice 

In this section the magnetic transitions of the Hubbard 
model on a simple-cubic (SC) lattice in three dimensions 
will be examined. The magnetic phase diagram of the 
Hubbard model has been studied with a multitude of 
methods for various lattices. These include perturbation 
theory in the interaction C/, 37 direct QMC, 38 ' 39 diagram- 
matic approaches, 40,41 DMFT, 12,42-46 and cluster exten- 
sions thereof. 47 

The reason we redo such an analysis here, is to in- 
vestigate possible shortcomings of the present approach 
by comparing its results to the literature. The impurity 
solvers based on the hybridization expansion are known 
to have some limitations in the Fermi liquid regime at too 
low temperatures. 22 ' 48 Additionally, the RPA-like decou- 
pling of the Bethe-Salpeter equations introduces a further 
approximation as envisaged above. As such decoupling 
schemes are known to overestimate the tendency toward 
phase transitions, we investigate the quality of this ap- 
proximation. 

The single-particle dispersion for the SC lattice with 
nearest-neighbor hopping in three dimensions is given by 



tk = — 2i[cos(fc x ) + cos(fc y ) + cos(fc 2 )] . 



(25) 



The half bandwidth is W = 6t and will be used as the 
unit of energy in the following. 

Figure 1(a) displays the temperature dependent in- 
verse static susceptibility for a half-filled SC lattice in 
three dimensions and various wave vectors q. Xmag(q, v = 
0) _1 depends almost linearly on T over the whole range 
of temperatures, which is expected from the mean-field 
character of DMFT. 12 ' 49,50 The wave- vector dependence 
of Xmag(?, v — 0) 1 is shown in Fig. 1(b) for various tem- 
peratures. The Q = 7r(l, 1, 1) T component of the suscep- 
tibility is always largest, indicating the strong tendency 
toward AFM. The dashed lines are fits to an expansion 
of the inverse susceptibility up to second order in the 
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FIG. 1. (Color online) Inverse static susceptibility (a) as func- 
tion of temperature for various values of the wave vector q and 
(b) as function of q for three different temperatures. Both 
plots are calculated for the SC lattice with nearest neighbor 
hopping only (t' = 0), U/W — 1.5, and at half-filling n = 1. 



hopping, 



Xfit(q,v = 0) 



1 



a{T) + b{TW) q 



(26) 



wher 



ere (i 2 ), = j^Y.^h.+^h.- For tne simple-, body- 
centered, and face-centered cubic lattice with nearest 
and next-nearest neighbor hopping, t and t' respectively, 
this elementary two-particle dispersion is just given by 
the negative single-particle dispersion, where all hopping- 
matrix elements are replaced with their squares, 



l\t^t 2 ,t> 



(27) 



(The sign stems from the definition of the single-particle 
dispersion with an overall negative sign, see Eq. (25.)) 
This form approximates the (^-dependence of the actual 
susceptibility very well, as it is visible in Fig. 1(b). 

This result provides a justification of the RPA-like de- 
coupling of frequency sums in the Bethe-Salpeter equa- 
tions. The reasoning is as follows: Comparing the q- 
dependency of the approximate form Eq. (26) to Eq. (17), 
we can infer that the particle-hole propagator - as the 
only g-dependent quantity entering the susceptibility - 



must be of the form 
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Po 



l-a(t 2 ) ( 



= P0 EW 



(28) 



n=0 



The geometric series, used in the last equality, can be 
interpreted as the result of an RPA-likc decoupling of 
the frequency sums appearing in the expansion of the 
exact particle- hole propagator, see Eq. (B8). Given the 
quality of the fit, we are led to the conclusion that the 
RPA-like decoupling works very well for the particle-hole 
propagator. It therefore seems reasonable to assume that 
it is also a good approximation in the Bethe-Salpeter 
equations. 

The physical reason lies in the nature of antiferro- 
magnetic correlations. They favor neighboring spins to 
point in opposite directions, which does not require a 
precise adjustment in the time-domain for a propagat- 
ing particle-hole pair. Therefore, decoupling in frequency 
space cannot induce a qualitative error (in contrast to the 
ferromagnetic case, see below). We always found good 
agreement between the form of Eq. (26) and the calcu- 
lated susceptibility, whenever the tendency toward AFM 
was dominant. 

For a larger Coulomb interaction of U = 3.3W and 
at finite doping, 8 = 3% (n = 0.97), the AFM suscep- 
tibility is largest only at high temperatures. This can 
be observed in Fig. 2(a) where x(9;0) -1 decreases lin- 
early with temperature only for T/W > 0.2 and changes 
qualitatively at low T. There, the FM component be- 
comes dominant. The wave-vector dependency as shown 
in Fig. 2(b) is well approximated by Eq. (26) for the 
AFM-dominated region at high T, but below T < 0.12W 
the fit does not work well and strong deviations occur. 

The AFM transition temperature Tn (Neel tempera- 
ture) is shown in Fig. 3 as function of the Coulomb repul- 
sion for the half-filled (8 = 0) model. Results from other 
studies arc also shown for comparison. The present result 
should be close to the data of Jarrcll (Ref. 12) or Ulmke 
et al. (Rcf. 43) as these results were also obtained with 
DMFT, but for a Gaussian (Jarrell) and semi-circular 
(Ulmke et al.) free DOS. For U < W the agreement is 
quite satisfactory, while for larger U the Neel tempera- 
ture from the present approach is substantially smaller. 
It is, however, roughly in agreement with data from QMC 
simulations of the three-dimensional model (Staudt et al, 
Ref. 38) but this may be viewed as a coincidence. 

The exchange interaction decreases as 1/U in leading 
perturbative order for large U. Moreover, it is renor- 
malized by small quasiparticle weights (Z-factors) when 
quasiparticlc bands emerge. So a decrease of Tn at larger 
U is to be expected, and its steepness depends on the size 
of the Z-factors, 



Z = 



1 



1 



dRcE 



(0) 



(29) 



While the ENCA as our impurity solver gives the cor- 
rect order of magnitude and parameter dependence for 
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FIG. 2. (Color online) (a) Inverse static susceptibility as func- 
tion of temperature for various values of the wave vector q. (b) 
1/Xma g (g, 0) as function of q for the temperatures (top to bot- 
tom): T/W = 0.26, 0.17, 0.14, 0.12, 0.065, and 0.021. Both 
plots are calculated for the SC lattice with nearest neighbor 
hopping only (t' = 0), U/W = 3.3, and n = 0.97. 



the low-energy scale of the SIAM, it underestimates its 
actual value. 22 Within the self-consistent treatment of 
DMFT this tendency is retained 52 as it can be observed 
in Fig. 4 where Z is shown as function of U. For compar- 
ison, the result from DMFT calculations with the NRG 53 
as impurity solver for the same parameter values is also 
shown. While the ENCA quasiparticle weight is smaller 
than the one extracted from the NRG, both follow the 
same trend. The inset displays the ratio of both and 
reveals the ENCA to yield a low-energy scale which is 
roughly a factor of 3 too small. One may argue, that the 
ENCA-calculation is performed for an effectively larger 
value of the interaction U. Therefore, the reduction of 
the critical temperature compared to the other DMFT 
calculations as observed in Fig. 3 is a result of the too 
small low-energy scale of the ENCA. 

Figure 5 displays the regions in the 8 — U phase-space 
where magnetic phase transitions are observed. The 
AFM region extends from half-filling 8 = to finite dop- 
ing for Coulomb interactions U < 2.3W (U/(W + U) < 
0.7). The investigation of regions in phase space with 
larger doping or smaller Coulomb repulsion was not pos- 
sible, as there the pathology of the ENCA 22 ' 48 prevented 
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FIG. 3. (Color online) The Neel temperature of the half- 
filled Hubbard model for the three dimensional SC lattice 
as function of U (red line with points). For comparison, re- 
sults from other studies are included as well: DMFT+QMC 12 
(blue crosses, labeled Jarrell), DMFT+QMC 43 (green stars, 
labeled Ulmke et al.), QMC 38 (brown squares, labeled Staudt 
et al.), high temperature expansion 51 (dash-dotted line, la- 
beled 3.83t 2 /U), and mean-field (dashed-line, labeled 6t 2 /U). 
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FIG. 4. (Color online) The quasiparticle weight for the half- 
filled Hubbard model as function of the Coulomb repulsion U 
extracted from the DMFT with ENCA and NRG as impurity 
solvers. The inset shows the ratio of both. 



FIG. 5. (Color online) Magnetic phase diagram in the S-U 
plane. Each point indicates a magnetic phase transition. (The 
regions without points were unaccessible with the ENCA as 
impurity solver). The inset shows the transition temperature 
as function of U for a fixed doping of 5 = 1 — n — 0.03. In 
both graphs, the red crosses indicate AFM, the blue circles 
IC, and the green squares FM transitions. 
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FIG. 6. (Color online) Inverse static susceptibility as func- 
tion of wave vector q for various values of temperature (top 
to bottom): T/W = 0.17, 0.12, 0.093, 0.056, 0.033, 0.016, and 
0.012. The curves are calculated for the SC lattice with near- 
est neighbor hopping only (t' = 0), U/W = 3, and n = 0.97. 



the access of low-enough temperatures. 

In accord with earlier studies 15 ' 37,45 ' 54 we also find in- 
commensurate (IC) spin-density wave transitions away 
from half filling at the edge of the AFM region. As a func- 
tion of U (see inset) or doping (not shown) , the transition 
temperature follows the trend foreshadowed by the Neel 
temperature, but the ordering vector shifts away from 
the AFM Q = 7r(l, 1, 1) T . 45 

This can be observed in Fig 6, where the inverse mag- 
netic susceptibility is shown as function of q for various 
temperatures. There the transition vector probably lies 



in the region around q = (1, 1, 0) T ir but a precise location 
is not possible since we did not scan a sufficient part of 
the Brillouin zone but only along selected axes. 

For finite doping, the nesting property of the Fermi 
surface is lost and incommensurate ordering is expected 
with increasing doping. 45 At a fixed doping, as presented 
in the inset of Fig. 5, the transition is induced by in- 
creasing the interaction. We attribute this to the tem- 
perature dependent formation of quasiparticlcs, whose 
characteristic energy scale is larger at lower U. There, 
the transition takes place at higher T where the details 
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of the Fermi surface are washed out and the nesting is 
approximately better fulfilled which leads to a AFM tran- 
sition. Increasing U decreases the transition temperature 
and the details of the Fermi surface become increasingly, 
making the system more sensitive to the lack of nesting 
and leading to the IC transition. 

As depicted in the inset of Fig. 5, the AFM and FM 
phases seem to be separated by a quantum critical point 
(QCP), where the transition temperature vanishes and 
the systems goes from AFM to FM ground states at 
T = 0. This QCP could in principle be physical as was 
speculated in the literature 55 but this question can not 
be addressed within the present study. The ENCA as im- 
purity solver does not allow the investigation of very low 
temperatures and consequently whether or not the tran- 
sition temperature vanishes cannot be decided. As FM 
is overestimated within the present approach (see below) 
it could well be shifted to larger values of U in more ac- 
curate treatments and thus removing the apparent QCP. 

At very large U > 3W we observe FM at finite doping, 
which is the region where it is in principle expected and 
observed for various lattices. 49 ' 56 However, for the un- 
frustratcd {t 1 = 0) SC lattice studied here, the transition 
temperatures Tq are too large (see inset). 

We attribute this overestimation of the tendency to- 
ward FM to the RPA-likc decoupling. Time-coherent 
propagation of particle-hole pairs over long distances is 
known to play a major role for FM 35 (ordering vector 
Q = 0). FM correlations favor equal spin-alignment, and 
a transfer of electrons between sites does require a time- 
correlation with accompanying holes due to the Pauli 
principle. But such correlations in the time domain are 
not conserved by the RPA-frequency decoupling at inter- 
action vertices. This conclusion is supported by the fact, 
that the FM transition temperatures as shown in the in- 
set are roughly in accord with the Stoner-like criterion, 

p(Q,Tc)U = 1, where p(0, T) is the fully interacting and 
temperature dependent DOS at the Fermi level. 

We have now established that the present approach 
yields reliable results whenever short-ranged particle-hole 
excitations dominate the magnetic response. In and close 
to the AFM regime, the transition temperatures are even 
reproduced quantitatively, apart from a reduction of Tn 
due to the too small low-energy scale produced by the 
ENCA. For very large Coulomb repulsions, the tendency 
towards FM is found to be overestimated which is at- 
tributed to the RPA-like decoupling of frequency sums 
in the Bethe-Salpeter equations. 



B. Body-centered cubic lattice 

In this section we focus on the three-dimensional body- 
centered cubic (BCC) lattice with nearest and next- 
nearest neighbor hopping, t and t', respectively. The 



dispersion is given by 

it = — 8£cos(fc x ) cos(ky) cos(fc z ) (30) 
- 2i'(cos(2Jfe x ) + cos(2fe y ) + cos(2fc 2 )), 

where the half bandwidth is W — 8t as long as \t'\ < 
which is always the case in the following. As for the SC 
lattice, the BCC lattice is bipartite and exhibits perfect 
nesting for vanishing next-nearest neighbor hopping, t' = 
0. The AFM nesting vectors Q satisfy tk+Q = — £fc and 
are of the type Q = tt(1, 0, 0) T and Q = tt(1, 1, 1) T . 

The perfect nesting property is illustrated in Fig. 7, 
where the Fermi surface of the nonintcracting system at 
half-filling is shown for two values of the next-nearest 
neighbor hopping, t' = and t' = —0.2t. The flat and 
parallel sections of the Fermi surface are clearly visible 
for t' = in panel (a). For finite t' (see panel (b)) the 
Fermi surface acquires a substantial curvature so that the 
nesting vector no longer connects large parts of the Fermi 
surface. 

The nonintcracting DOS as shown in Fig. 8 displays 
the characteristic van-Hove singularities due to the ex- 
trema and saddle points of the dispersion relation tk in 
the Brillouin zone. Most prominent is the logarithmic 
divergence near the Fermi level due to the maxima at 
k = -|(1, 0, 0) T . Such a van-Hove singularity can induce 
profound changes in the low temperature, low-energy 
Fermi liquid properties, even in DMFT. But here we 
do not study this possibility but focus on magnetic tran- 
sitions which generally occur at higher temperatures. 

The reason we are considering the BCC lattice instead 
of the more common SC lattice is found in the more pro- 
nounced influence of the frustrating next-nearest neigh- 
bor hopping t' in the former lattice. The shape of the 
Fermi surface changes much more strongly with increas- 
ing t' for the BCC lattice as, for example, in the SC 
lattice. In particular, results pointing to a competition 
between local-moment and quasiparticle magnetism are 
found to be more pronounced for the BCC lattice than 
for the SC case, even though we also observed them in 
the latter case as well. 

A DMFT-spcctral function for the BCC lattice with fi- 
nite next-nearest neighbor hopping is shown in Fig. 9. As 
a consequence of the rather large Coulomb interaction the 
lower and upper Hubbard band are clearly visible. Upon 
lowering the temperature a quasiparticle peak emerges at 
the Fermi level indicating the formation of a correlated 
metallic state. The inset shows the imaginary part of the 
self-energy. The quadratic minimum at the Fermi level 
reveals this state to be a Fermi liquid as it is expected 
within DMFT. 11:57,58 The momentum resolved spectral 
function p(k, oS) shown in panel (b) exhibits the dispersive 
low-energy quasiparticle band around the Fermi level. 

The formation of this low-temperature Fermi liquid is 
associated with a screening of local magnetic moments 
due to Kondo-correlations. Figure 10(a) displays the 
screened local moment, 

Mcff = T X l Z s (T), (31) 
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FIG. 7. (Color online) Noninteracting Fermi surface for the 
BCC lattice for (a) t' = and (b) t' = -Q.2t at half-filling. 
(In the front quadrant some parts are omitted in order to get 
a better view.) 

as function of temperature for three different values of the 
Coulomb interactions. Due to the Coulomb interaction 
the moments at high temperatures are larger than the 
induced moment /j,^ s u=0 = | of the noninteracting band 
electrons (where empty and doubly occupied local states 
fully contribute), and increases with U toward the value 
Meff spin = I °f a f ree s P m - At temperatures of the order 
of the characteristic low-energy scale T* = WZ, where 
Z is the quasiparticle weight of Eq. (29), the moments 
decrease due to the buildup of Kondo-corrclations. 
In contrast, the unscreened moment 

Mbaxe - \ [<"t> + {H) - 2("t"i>] (32) 

which is essentially determined by the double occupancy 
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FIG. 8. (Color online) Noninteracting DOS for the BCC lat- 
tice at half-filling and three different next-nearest neighbor 
hoppings, t'/t = 0, -0.1, and -0.2. 



is almost independent of temperature, see Fig. 10(b). 
Upon lowering the temperature it slightly increases due 
to the decrease in the thermally induced double occu- 
pancy. The formation of the Fermi liquid quasiparticle 
band at low temperatures then leads to a slight increase 
of the double occupancy 40,59 and consequently the mo- 
ment decreases again. 

AFM should prevail for not too large next-nearest 
neighbor hopping t' due to the near-nesting property of 
the Fermi surface. Figure 11(a) displays the inverse mag- 
netic susceptibility for U/W = 1.25 and t'/t = —0.15 
as function of the wave vector along a path through 
the Brillouin zone. As in the previous section for the 
SC lattice, the AFM components (Q = 7r(l,0,0) T and 
Q = 7r(l,l,l) T ) are enhanced. Additionally, fits with 
the approximate function of Eq. (26) agree very well 
with the susceptibility, supporting the view developed 
in the previous section. The temperature dependence of 
the susceptibility is depicted in Fig. 11(b). The AFM 
component exhibits a Curie- Weiss behavior T ^ T ^ with a 
divergence at T^/W « 0.035, which supports an inter- 
pretation in terms of local-moment magnetism. 

However, this perspective becomes less obvious when 
the constituents are analyzed in detail. The effective lo- 
cal moment which enters Eq. (17a) is temperature de- 
pendent due to screening, see Fig. 10. Additionally, the 
static exchange coupling J(q, v = 0) also displays a tem- 
perature dependency. This can be observed in Fig. 12(a), 
where the AFM coupling J(Q,0) slightly increases to- 
ward lower temperatures. The increase in J(Q, 0) is com- 
pensated by the reduced moment and together both yield 
the observed Curie- Weiss behavior characteristic of local- 
moment magnetism. 

Figure 12(b) shows the same data but organized in 
terms of the itinerant picture of magnetism. The "nonin- 
teracting" susceptibility, i.e. the particle-hole propagator 
calculated without explicit two-particle interactions but 
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temperature for t'/t = —0.15 and three different values of the 
Coulomb interaction. 



FIG. 9. (Color online) (a) Local spectral function p(uj) ob- 
tained with the DMFT for the BCC lattice for three different 
temperatures. The curves are calculated for half-filling with a 
Coulomb interaction U/W = 1.375 and a finite next-nearest 
neighbor hopping t'/t = —0.2 is used. The inset shows the 
(negative) imaginary part of the self-energy -ImE(cj + iS) for 
the lowest temperature T/W = 0.01. (b) The momentum de- 
pendent lattice spectral function as function of k along high 
symmetry directions in the Brillouin zone for the same pa- 
rameters as in (a) and T/W = 0.01. 



with the full single-particle Green functions, increases 
with decreasing temperature for all wave vectors. This is 
a consequence of the accumulation of spectral weight near 
the Fermi level, cf. Fig. 9(a). The effective magnetic two- 
particle vertex T miig (i> = 0) is of the order but slightly 
larger then the bare interaction U (see inset). It slightly 
decreases for low temperatures indicating the quasipar- 
ticle excitations to experience somewhat less scattering. 

Thus, the buildup of Kondo-like resonances in the lat- 
tice has two opposing effects, which nearly compensate 
each other here: The screening of local moments and the 
accumulation of additional quasiparticle weight near the 
Fermi level. 

The balance between itinerant and localized contribu- 



tions to the susceptibility is rather fragile. Already for 
moderate changes of parameter values the Curie- Weiss 
behavior of the susceptibility is lost. This is the case 
for a slightly larger next-nearest neighbor hopping of 
t'/t = -0.17 and U/W = 0.875 as shown in Fig. 13(a). 
AFM is still dominant but no magnetic transition occurs. 
All components of the susceptibility exhibit a maximum 
(minimum in inverse susceptibility) and decrease toward 
lower temperatures. The decrease occurs at a tempera- 
ture of the order of the low-energy scale indicating that 
the coherent quasiparticles do not favor any tendency to- 
ward magnetic order. 

This maximum in the susceptibility opens up the pos- 
sibility of an unusual phase-transition scenario. Increas- 
ing the Coulomb repulsion enhances the AFM suscep- 
tibility as can be seen in Fig. 13(b). For large enough 
U, the minimum in l/x m ag(Qj0) even reaches zero in- 
dicating a divergent susceptibility at a finite tempera- 
ture (Tn/W w 0.016 for the case shown in the graph 
with a Coulomb repulsion between U/W — 1.025 and 
U/W = 1.05). 

Fixing the temperature at some value and increasing 
U can lead to the peculiar situation that the system stays 
paramagnetic even though at a higher temperature it 
would undergo an AFM phase-transition. In that regime, 
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FIG. 11. (Color online) (a) The inverse magnetic susceptibil- 
ity for different temperatures as function of the wave vector 
along with fits of the form a + b{t 2 ) q , see Eq. (26). (b) The 
inverse magnetic susceptibility for different wave vectors as 
function of temperature. The curves are calculated for half- 
filling with a Coulomb interaction U/W = 1.25 and a next- 
nearest neighbor hopping of t' /t — —0.15. 



FIG. 12. (Color online) (a) The effective static exchange cou- 
pling J(q, v = 0) of Eq. (18) for different wave vectors as 
function of temperature for U/W = 1.25 and t'/t — —0.15. 
(b) The inverse negative static particle-hole propagator for 
different wave vectors as function of temperature. The inset 
displays the effective local vertex of Eq. (22). 



there exists a finite temperature interval Tn.i <T< Tn ; 2 
in which AFM order is expected to occur, while outside 
this interval the system stays paramagnetic. Unusually, 
for T < Tn.i the system can be driven into the AFM 
ordered phase by heating it. 

It should be noted, that similar re-entrant behavior 
was already found for the AFM transition in the Hubbard 
model within a weak-coupling treatment for an infinitc- 
and two- dimensional lattice, 60,61 as well as within an 
equation of motion technique. 62 In these cases, the frus- 
tration was brought about by finite hole-doping, instead 
of a finite next-nearest neighbor hopping implemented 
here. The extended Hubbard model with an additional 
nearest-neighbor Coulomb interaction also displays a re- 
entrant behavior for the charge-ordering transition. 63 

Another peculiarity feature is connected to the critical 
exponent. The theory we employed has mean- field char- 
acter regarding spatial correlations. Thus, we expect a 
mean-field critical exponent 7 = 1 for the nonlocal mag- 
netic susceptibility, x{Q, T ) ~ ( T ~ T c )~~< ■ This expo- 
nent is indeed always found, except for that value of U 
where the minimum in touches the zero-axis. There, 



the exponent is is larger and of the order 7 = 2 as this 
minimum can be fit with a parabola. 

The re-entrant behavior is a result of the competition 
between the ordering of local moments and the forma- 
tion of low-temperature quasiparticles, i.e. the compe- 
tition between local moment and itinerant quasiparticle 
magnetism. The exact shape of the Fermi surface has 
a strong influence on the quasiparticle scattering which 
determines the particle-hole propagator. At low tempera- 
tures, the quasiparticles are fully formed and the frustra- 
tion induced by t' leads to a suppression of the particle- 
hole propagator, essentially eliminating magnetism (sim- 
ilar to what happens at weak coupling). 

This can be observed in Fig. 14(a), where the in- 
verse particle-hole propagator shows a saturation and 
even a slight upturn toward low temperatures (compare 
Fig. 12(b) where this is not the case), even though the 
system is metallic with large spectral weight at the Fermi 
level, cf. Fig. 9(a). The effective vertex (see inset) is 
not sensitive to the Fermi surface and consequently dis- 
plays a qualitatively similar behavior to the case shown 
in Fig. 12(b). 
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FIG. 13. (Color online) (a) The inverse magnetic susceptibil- 
ity for different wave vectors as function of temperature for 
U/W = 0.875. (b) The inverse AFM component of the sus- 
ceptibility for small temperatures and various values of U. In 
both plots is t'/t = -0.17. 



A local moment perspective can be employed above 
T/W « 0.02 where the effective exchange coupling J 
is almost temperature independent (see Fig. 14(b)) and 
very similar to the one shown in Fig. 12(a). Below this 
temperature, however, coherent quasiparticles emerge 
and the coupling changes dramatically. The AFM com- 
ponent even changes sign and becomes ferromagnetic. 

Whether the local moments order at higher tempera- 
ture is determined by their size and the geometry of the 
lattice. Both aspects are much less sensitive to the frus- 
tration generated by a finite t' than the magnetic quasi- 
particle scattering. The value of the screened effective 
moments is not influenced much by t' as can be seen in 
Fig. 15(a). /^ ff is shown for a fixed temperature as func- 
tion of t' for two different values of U. The decrease with 
increasing t' is barely visible confirming the screening to 
be independent of the details of the Fermi surface or the 
shape of the spectral functions (as long as the system is 
a metal, of course). 

In order to get an estimate for the t' dependence of the 
Neel temperature from the local moment perspective, it 
is instructive to start from the atomic limit. The effec- 
tive exchange coupling can be estimated in second order 



FIG. 14. (Color online) (a) The inverse negative static 
particle-hole propagator for different wave vectors as func- 
tion of temperature. The inset displays the effective local 
vertex of Eq. (22). (b) The effective static exchange coupling 
J{q,v = 0) of Eq. (18) for different wave vectors as function 
of temperature, for In both plots U/W = 1 and t'/t = —0.17 
is used. 



perturbation theory in the hopping, 
8t 2 

Jeffig) = — — cos(q x ) cos(q y ) cos(q z ) 



(33) 



2T 
IT 



cos(2(7 x ) + cos(2g y ) + cos(2g 2 ) 



which yields the effective AFM coupling for q = Q. 
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3^ 
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This coupling directly determines the AFM transition 
temperature as it was discusses around Eq. (21). The 
Hubbard model for the parameter values under consider- 
ation is far from the atomic limit as it is a metal and 
itinerant electronic excitations are present. However, 
for the local moment picture to be applicable the Neel 
temperature should follow the quadratic decrease with t' 
of Eq. (34), but renormalized pre- factors would be ex- 
pected. Figure 15(b) displays Tn as function of t' for 
two different values of U . The transition temperature 
indeed decreases in accordance with that expectation for 
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FIG. 15. (Color online) (a) The effective local moment as 
function of the next-nearest neighbor hopping t' for a fixed 
temperature T/W = 0.046 and two values of the Coulomb 
repulsion, (b) The Neel temperature as function of t' normal- 
ized to its value for t' = for two values of U. Also shown are 
fits with a quadratic decrease, where the fit parameters are 
a = 24 and a = 10.95 for U/W = 0.875 and U/W = 1.375, 
respectively. 




U/W 



FIG. 16. (Color online) The Neel temperature of the half- 
filled Hubbard model on BCC lattice as function of U for 
various t . 



AFM can occur. Similar behavior is found for the SC lat- 
tice with finite next-nearest neighbor hopping i'. 44 The 
perturbation expansion in U shows that for vanishing 
next-nearest neighbor hopping, t' = 0, AFM extends to 
[7 = 0, i.e. Uc = 0, due to perfect nesting (see, for exam- 
ple, Ref. 1, Chap. 7). 

For large U AFM is found to be very robust with re- 
spect to t' and Tn is only reduced. The dependence of the 
Neel temperature on t' is in agreement with the quadratic 
decrease expected from a local-moment picture, and the 
quadratic coefficient of the t'-term decreases with increas- 
ing U, cf. Fig 15(b). In the crossover regime at interme- 
diate U/W the re-entrant behavior is observed. For a 
given U, there exist two critical temperatures, Tn,2 and 
Tn,i, between which the system exhibits AFM, while be- 
low Tn,i and above Tn,2 paramagnetism prevails. 



small to moderate values of \t'/t\. The rather large val- 
ues for the fit parameter a compared to 3/4 of Eq. (34) 
indicate strong rcnormalizations due to the formation of 
heavy quasiparticles. At too large t', the Neel temper- 
ature comes in regions where the forming quasiparticles 
begin to dominate and magnetism is suppressed quite 
abruptly. For smaller values of U, the local moments are 
generally less pronounced and consequently this suppres- 
sion occurs faster. 

The Neel temperature does not go to zero continuously, 
but rather remains finite at the maximal t 1 where a tran- 
sition occurs. In the vicinity of this value the re-entrant 
behavior can be observed, and the system tends to ex- 
hibit AFM only in a finite temperature interval. 

The Neel temperature as function of the Coulomb re- 
pulsion U for various values of the next-nearest neighbor 
hopping t' is shown in Fig. 16. With increasing t' the 
transition temperature is generally reduced especially in 
the weak coupling regime where AFM is completely elim- 
inated. For finite t' , there exists a nonzero critical value 
for the Coulomb repulsion, Uc > 0, only above which 



IV. CONCLUSION AND OUTLOOK 

We have derived an universal and physically intuitive 
form of the magnetic susceptibility of strongly correlated 
electron systems within DMFT using an additional RPA- 
like decoupling in the Bethe-Salpeter equations. It was 
thus possible to work completely in the real-frequency do- 
main and to write this final expression in terms of func- 
tions of real variables, which all are obtained from an 
effective impurity model. The static and dynamic mag- 
netic response could readily be calculated from there. 

We analyzed the validity of this form and found 
the RPA-like decoupling to be quantitatively correct 
whenever short-ranged particle-hole excitations domi- 
nate. This is the case in and close to the antiferromag- 
netic region in phase space of the Hubbard model where 
the Neel temperature was correctly reproduced (modulo 
shortcomings of the ENCA impurity solver). 

The tendency toward fcrromagnctism, on the other 
hand, is overestimated by the approach as a direct con- 
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sequence of the decoupling. However, this shortcoming 
could be circumvented in the future by employing a dif- 
ferent decoupling strategy, where averages are taken over 
products of a particle-hole propagator and a local vertex. 
Leading order coherence effects during a two-particle in- 
teraction should then be retained. 

The applicability of the two archetypical pictures of 
magnetism was then investigated in case of the Hubbard 
model on a three-dimensional body-centered cubic lattice 
with next-nearest neighbor hopping. For small interac- 
tions strength the antiferromagnetism can be understood 
in terms of itinerant quasiparticlcs. The magnetic re- 
sponse is determined by particle-hole excitations across 
the Fermi surface. With increasing next-nearest neigh- 
bor hopping t' , the perfect-nesting property of the Fermi 
surface is lost rapidly. As a consequence, AFM is almost 
entirely removed from the weak coupling phase diagram. 

At large values of the local Coulomb interaction U the 
local-moment picture is more appropriate, as the anti- 
ferromagnetism is rather robust against geometric frus- 
tration induced by t' . The expansion around the atomic 
limit provides an suggestive form for the t'-dependence 
of the Neel temperature which is in agreement with the 
DMFT results. 

At intermediate U both types of magnetism compete. 
The many-body effects leading to the formation of the 
quasiparticlc bands and the screening of local moments 
are strongly temperature dependent. This produces an 
unusual re-entrant behavior. Local moments still sizable 
at elevated temperatures order antiferromagnetically at a 
characteristic temperature Tn,2- But at even lower tem- 
peratures the coherent quasiparticles are dominating and 
due to the lack of nesting of the Fermi surface, magnetism 
is suppressed below Tn,i- 

Beyond the scope of magnetism, the interplay of lo- 
calized and itinerant effective pictures is of fundamental 
interest. For example, the explanation of the re-entrant 
Mott-transition found in the highly frustrated organic 
compound k-(ET)2Cu[N(CN)2]C1 7 is in direct analogy 
to the mechanism presented here. The temperature de- 
pendency of the effective picture and the re-entrant be- 
havior is generated by a competition between frustrated 
magnetic moments and a low temperature Fermi liquid. 

It is an interesting prospect and needs further investi- 
gations whether similar temperature dependencies could 
play a role in other strongly correlated systems. In 
the cuprate superconductors, 3 for example, the single- 
particle excitations are believed to be of a more itiner- 
ant character along the nodal direction while of more 
localized nature along the anti-nodal. 64 The frustration 
brought about by the next-nearest neighbor hopping t 1 
could play a decisive role, at least for some compounds 
and parameter values. 
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Appendix A: DMFT for single particle Green 
functions 



In this section we sketch the DMFT self-consistency 
which will turn out to be useful for the formulation of 
the two-particle properties in appendix B. The Dyson- 
equation can be written in terms of the effective local 
single-particle cumulant Green function 



G(z) 



1 



(Al) 



This quantity represent dressed atoms 65 between which 
single-particle transfers are possible. This view is es- 
tablished within the linked cluster expansion around the 
atomic limit utilizing local cumulants 22 ' 27,65 ' 66 and yields 
the Dyson equation 



G(k, z) = G{z) + G{z) t k G(k, z). 



(A2) 



Solving this equation for the Green function just gives 
Eq. (4). 

The local Green function is obtained by a Fourier 
transform 



1 



G(z) = G(z) + G(z)jjY,tkG(k,z) (A3) 

k 



G{z) + G(z) T{z)G{z). 



(A4) 



It was used that = 0, which holds for all inversion 
symmetric lattices and we introduced the local scattering 
matrix 



(A5) 



A graphical representation of this equation is depicted 
in Fig. 17(a). T(z) incorporates all processes where an 
electron leaves the site under consideration, propagates 
through the entire lattice — including that local site — 
and then returns. These processes can be reorganized in 
terms of re- visits to the local site, 67 



T{z) = T{z) + T(z)G{z)T{z) + 

= T{z)+T(z) G(z)T(z) 
1 



T(z)-i-G(z) 



(A6) 
(A7) 

(A8) 
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FIG. 17. (a) Dyson equation for the local single-particle 
Green function in terms of the local cumulant G (grey shaded 
circles) and the local scattering matrix T(z). The double line 
with the arrow and the small circles at the end indicate a 
full Green functions G(k, z) and a dashed line an elementary 
hopping tfc. (b) The local scattering matrix T(z) in terms 
of the irreducible loops, i.e. the effective medium, F(z) as in 
Eq. (A8). The full scattering matrix is already graphically 
composed of two elementary hoppings and a full propagation, 
see Eq. (A5). The effective medium is represented in a sim- 
ilar way, only the propagation through the lattice without 
visits to the local site is indicated by a zig-zag line. In both 
quantities the elementary hoppings are included and, there- 
fore, do not appear in Eq. (A6) - (A8). The start and end 
points of all propagations are connected to indicate local i.e. 
fc-summed quantities and the small horizontal bars indicate 
that these lines are not accounted for in an analytical expres- 
sion. The small circles are drawn only to indicate a lattice 
site at which a propagation starts/ends and do not appear in 
any analytical expressions. In the present context these are 
unnecessary, but they are helpful for two-particle quantities 
which are discussed in the next section. 



The effective medium T(z) now has the clear interpreta- 
tion of being the irreducible loop for a propagation of an 
electron from a specific site through the lattice back to 
that site, without returning during the propagation, see 
Fig. 17(b). This makes it very obvious to interpreted T as 
an noninteracting effective medium and thus establishes 
the DMFT mapping of the lattice problem to an SI AM. 

Formally, this is obtained by inserting Eq. (A8) into 
Eq. (A4) which gives the local Green function 



G(z) 



G~ 



T(z) 



(A9) 



which is nothing but Eq. (6). 



Appendix B: Bethe-Salpeter equations and dynamic 
susceptibilities 

The physical susceptibility matrix depends on one ex- 
ternal bosonic Matsubara frequency iv n and wave vector 
q. It is obtained from a more general two-particle Green 



FIG. 18. Graphical representation of the susceptibility matrix 

Xa,b;b,a{q, W n \k x , fc 2 \w\ , IU2) . 



FIG. 19. Bethe-Salpeter equation for the susceptibility ma- 
trix. The shaded square represents the irreducible vertex II. 



function by summing over two internal fermionic frequen- 
cies and wave vectors (omitting the convergence factors 



BN 

iuji ,iaj 2 k-y ,k 2 
= 7? 51 X(9,^n|«Wi,iW2). 



(Bl) 
(B2) 



Here, iu>i and iw2 are short-hand notations for the sets of 
fermionic Matsubara frequencies ico ni and iu> n2 i n i, n 2 £ 
Z). The notation is introduced to discriminate between 
external arguments (iu n , q) and internal momenta {k 11 
k 2 ) and frequencies {iw\, IW2) which are separated by 
vertical bars. In the last equation we absorbed the mo- 
mentum sum into the susceptibility as indicated by the 
missing internal momentum variables. 

We also introduced a matrix notation for two-particle 
quantities in orbital/spin space indicated by the straight 



double underline, i.e. A 



a,b:c,d 



{A}a 



b-c-d- 



The matrix 



multiplication is denned as 



{AR\a 



b:c,d 



ef 



A, 



f ,b:c,e^a,e;f,d' 



(B3) 



The graphical representation of the susceptibility ma- 
trix x(<Z> i v n 1^1,^2 ^2) is shown in Fig. 18. The 

two external pairwise contracted lines signal the specific 
choice of orbital matrix elements Xa,b;b,a, as well as fre- 
quency and momentum arguments. 

The Bethe-Salpeter equation for the susceptibility ma- 
trix is obtained by summing all particle-hole reducible di- 
agrams not accounted for in the vertex n. This is shown 

graphically in Fig. 19. 

Within DMFT the full particle-hole irreducible two- 
particle vertex II is assumed to be momentum indepen- 
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FIG. 21. Equation for the two-particle cumulant transfer 
propagator £ c . 



FIG. 20. (a) The particle-hole propagator expressed through 
Q and A GG . (b) The quantity A GG where the external lines 

are made explicitly for clarity. 



dent and momentum conservation is ignored at internal 
vertices. This allows all sums over internal momenta to 
be absorbed into effective quantities only depending on q. 
The particle-hole excitation of such an effective quantity 
thus always propagates between two sites only. Graphi- 
cally this is indicated in Fig. 19 by the endpoints of the 
particle-hole excitation always being right on top of each 
other. 

Analytically, the Bethe-Salpeter equation depicted in 
Fig. 19 is 



(B4) 



+ £(<?! iv n I w%) n aJnp (iw n I Vj3\ y iu3%) R(q, iv n \ 

The particle-hole propagator P_ is constructed from the 
single-particle Green function 



E(q, iv n | iui) = — GG( q, iv n \k\iwi) 



(B5) 



The matrix notation on the right-hand side indicates a 
two-particle matrix build from one-particle quantities by 
means of a tensor product 



GG (q, ii/ n \k\iu3i) 



(B6) 



Gaciki G b d(h + q, iuii + w n ) 
= G(k,iuJi)®G{k + q,ioJi + iv n ) (B7) 

In the last line we introduced the double under-wave as 
indicator for a one-particle matrix with two indices only. 

Inserting the Dyson equation (A2) and using 5Z fc ijt = 
0, the particle-hole propagator can be separated into a 
momentum-dependent and momentum-independent part 

P(£,ii/„ | iwi) (B8) 

= Miv n | iuJi) GG (iv n \ + Q(q,iu n | iwi), 



with the momentum dependency fully accounted for by 
the function 



Q(q, iv n | iui) 



(B9) 



iv n \k\iu)\) tt(a | k) GG iiVn \ iuii ) 



and 



Miv n | iu>i) =1+1 ® \_G(iuji+iv n ) T{iu)x+iu n )\ 



+ [T(iwi) G(iui)] ®1. 



(BIO) 



This decomposition is show in Fig. 20. The func- 
tion A GG is represented graphically as two circles with 

crosses and incorporates processes, where only one elec- 
tron leaves a site and propagates through the lattice, 
while the other remains at this site. 

The negative signs in front of the particle-hole propa- 
gators in Eq. (B4) stem from the closed Fcrmion loop, i.e. 
from the permutations necessary to obtain the ordering 
of the creation- and annihilation operators. Such signs 
are viewed as part of the diagrammatic rules and thus do 
not explicitly occur in the figures. 

In the Bethe-Salpeter equation (B4) the amputated 
version of the local vertex is used which is obtained by 
factorizing the one-particle contributions of the local site, 



(Bll) 



GG {iv n | iwa) W mp (iv n | iwi, iwi) GG{ iv n \ iui). 

This is done in order to avoid over-counting of contribu- 
tions from the start- and endpoint of the propagation. 
Graphically, amputation is indicated by the vertical bars 
around the arrows. 

In Eq. (B4) an additional quantity £ c had to be in- 
troduced in order to avoid two local irreducible vertices 
II occuring at the same site right after each other. This 
two-particle cumulant transfer-propagator describes cor- 
related particle-hole propagations starting with two ele- 
mentary transfers. A separate Bethe-Salpeter equation 
can be formulated 

£[ c (q, iv n | iwi,iu)2) = (B12) 
E(q, iv n | iw 2 ) W xap (iu n \ iuji,iuj 2 ) Q(q, iv n \ iui)- 

^2^{q-,w n | iuj3,icj 2 ) n amp (w„ | mji,kj 3 ) Q{g,iv n | iu>i), 
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FIG. 22. Alternative Bethe-Salpeter equations for the lattice 
susceptibility expressed with the effective local two-particle 
cumulant Green functions. 



which is shown graphically in Fig. 21. The crossed box 
signals, that only terms with at least one local irreducible 
interaction vertex II involving all four external lines con- 
tributes. 

Within DMFT, the local irreducible two-particle ver- 
tex II can in principle be determined from the effec- 
tive impurity model. But it is more instructive to to 
re-arrange the Bcthc-Salpcter equation (B4) in terms of 
dressed local two-particle Green functions, analogous to 
the one-particle Dyson equation (A2). The lattice sus- 
ceptibility can then be expressed as 



(B13) 



GG {iv n | ilo 2 ) S amp (q, iv n \ iuj 3 , iu 2 ) 



where the amputated two-particle transfer propagator 
gamp j s introduced. This equation is depicted in Fig. 22. 

The (un-amputated) two-particle transfer propagator 

is closely related to the two-particle cumulant transfer 
propagator S c via 



iv n | iu>i,Kd 2 ) = -Q{g, iv n \ ioJi)Siu uiUa (B14) 
The dressed local two-particle Green function entering 



Eq. (B13) is 

& ,A (iv n | iwx,iu)2) = -QQiivn l^i) <W,to> 2 ( B15 ) 
+ Il{ivn | iwi, zwa) Miv n I iw\ ) , 

which is shown graphically in Fig. 23. It is essentially 
determined by the local irreducible vertex but additional 
pseudo-local corrections due to one-particle propagations 
do occur. The minus sign associated with the first 
term reflects the fact that this term originates from the 
particle-hole propagator. Since the closed fcrmion loop 
is no longer visible in the diagrammatic representation it 
is included explicitly in the graphical representation. 
The factor of A appearing at the cumulant vertex in 

Eq. (B15) is not valid by the rules of the original cumu- 
lant perturbation theory, as two local contributions are 

cfW^p = -O + Q^ST^ 



= -H + 





^ + ^ • : 

, I-M- l>l 



FIG. 23. The effective local two-particle cumulant Green 
functions Q 2 , £ 2 ' A , AG^A" 1 (see footnote 68), and AG 2 ' A 

(from top to bottom). 



multiplied at the very same lattice site. The reason it ap- 
pears here lies in the internal momentum sums already 
performed in the lattice susceptibility iv n \ iuii,iu2) 

(see Eq. (B2)) together with the approximation of a mo- 
mentum independent vertex. Separating pseudo-local 
and nonlocal contributions the momentum-independent 
part of the particle-hole propagator inevitably leads to 
the appearance of such terms which are originally not 
allowed. 

Equation (B13) is the two-particle analog of the Dyson 
equation (A2) with 1PI diagrams replaced by 2PI dia- 
grams, along with the resulting factors of A, and the 

complication of a second equation needed for the two- 
particle transfer propagator 

For the mapping onto the effective impurity model we 
construct the local two-particle susceptibility by iterating 
Eq. (B13) once and summing over q which gives 68 
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g 

= Miv n 1 10J2) GL' K {iv n | iu)\, iui 2 ) 
+ ^2 M}v n | «w 2 )^ 2 ' A (w„ | iu) 2 )A(ii> n I iLOi)~ l TT{iv n \ iw 3 , iuja) £P' A (iv n | 



(B16) 




FIG. 24. Bethe-Salpeter equation for the local two-particle 
Green function x loc - 




FIG. 25. The decomposition of the two-particle scattering 
matrix TT into irreducible propagations TT and repeated 
visits to the local site. 



This equation is represented in Fig. 24. The quantity 
TT is closely related to the g-summed two-particle trans- 
fer propagator and represents the two-particle scattering 
matrix, where both particles leave the local site and prop- 
agate through the lattice. Equation (B16) is the analog 
of the one-particle equation (A4) and establishes the con- 
nection between the effective local two-particle cumulant 
Green function Q?' A (or equivalently the local cumulant 

vertex |J_ C ) and the local physical Green function. 

In complete analogy to the one-particle case (see 
Eq. (A8)) the two-particle scattering matrix TT is build 

up from an irreducible scattering matrix TT . The latter 

describes the simultaneous propagation of two particles 
starting and ending at the same site without any inter- 
mediate returns to this specific site. This decomposition 



is shown graphical in Fig. 25 which reads analytically 
TT_{iv n | iu)i, iu 2 ) = ^ I iuJ ?) <Ww 2 ( B17 ) 

+ ^] TT {iv n | iuji, iuj 2 ) Cl ,A (w,i I iw3,MJ 4 ) x 

x TT (iv„ I iui, 1LU3) h_{iv n \ ioji). 

The structural difference of this equation to the one- 
particle case Eq. (A8) lies in the appearance of the factors 
A which account for the single-electron recurrence. 

At this point the identity 



TT_{iv n \iu> 3 ,iu> 2 )G_ ' (iv n \ iw-i, iwz) 



(B18) 



= A(w„ I iw 2 )TT(w„ I iuj 3 ,iuj 2 )^ oc (iv n I iui,iuj 3 ) 

1U13 

and the equivalent form with the order of the matrix 
products reversed (as well as the order of the factors 
A) can be proven by comparing Eqs. (B16) and (B17). 

These reflect the two possible ways to decompose the lo- 
cal Green function and the effective local cumulant in 
terms of reducible and irreducible loops. 

Utilizing the above identities in Eq. (B16), the scatter- 
ing matrix can be eliminated in favor of the irreducible 
scattering matrix and the local susceptibility, 



M,iv n I iuJ 2 )Q?' A (iis n I iwi,KJ 2 ) 

X loc (w„ |«w 3 ,iw 2 )Tr(w„ \ i0j4, «w 3 )x 



(B19) 



E 

2UJ3 ,20J4 



x K(iv n 1 1UJ4) G_ ' iiv n I iui, iuj^). 

All equations derived so far incorporate inner fre- 
quency sums and thus have the character of integral equa- 
tions. The local susceptibility x loc is determined from so- 
lution of the impurity model and then used as an input 

_ 2 A 

functions for Eq. (B19). From this, G ' is determined 
which then has to be used to solve the lattice equations, 
Eqs. (B13) and (B12). Thus, closed explicit expressions 
for the momentum dependent susceptibilities cannot be 
obtained, in contrast to the single-particle case. 
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In order to complete the mapping onto the effective 
impurity model we need to specify the processes to be 
incorporated into the effective medium. Therefore, it is 
beneficial to analyze the structure of this set of equations 
by switching to the Matsubara-matrix notation. Equa- 
tion (B19) can be solved in Matsubara-matrix space to 
give 



-loc 



(B20) 



This clearly demonstrates the resemblance to the one- 
particle equation (A9). The dressed local two-particle 
Green function is given by AG 2 ' A , where the one-particle 

recurrence problem leads to the unexpected factor of A, 
and TT is the irreducible effective two-particle medium. 

The effective local two-particle vertex can now be ex- 
pressed in terms of the local susceptibility function and 
used to obtain an explicit form for the lattice suscepti- 
bility in Matsubara-matrix space 



Xq,- 



(B21) 



As the effective impurity is embedded into an non- 
interacting medium within the DMFT, the irreducible 
two-particle medium TT can be specified in terms of the 
effective medium for the one-particle Green function, 

T~T(iv„ I iui-i , iuj?) = (B22) 

- (T(iui) ® T(iui + iv n )) Miv n | iwi) -1 <Li,w a - 

The negative sign reflects the fact that the effective 
medium is derived from the g-summed particle-hole prop- 
agator, where a minus sign has to be included due to the 
closed fcrmion-loop. The factor of A" 1 removes those 
processes where only one electron intermediately re- visits 
the site, which are present in the product of the two one- 
particle media T. This is necessary as TT characterizes 

the amplitude where both electrons must simultaneously 
leave the local site. 

Notice that in this approximation the scattering ma- 
trix TT is not just the product of the uncorrelated one- 
particle scattering matrices, i.e. 

TT_{iv n | iivi,iLo 2 ) ^ ^T( iuj i) ®T(iuji + iv n )6 UuUJ21 

(B23) 

as it could be suspected naively. The reason is that two- 
particle correlations are discarded only in the medium, 
but locally all interaction vertices are retained. The 
Matsubara-matrix space representation reveals the struc- 
ture of TT 



This very instructive form confirms the insight, that the 
uncorrelated one-particle loops get rcnormalized by in- 
teraction vertices whenever the particles visit the local 
lattice site. 

The explicit form (B22) can be used to express the 
irreducible scattering matrix as 



(B25) 



— \Miv n | Mi ) GG (iv n I iuii) 
where the local particle-hole propagator is 

P loc ( W „|z Wl ) = -^J^P(g,iVn\iui) (B26) 



= G 1oc (mj!) <g> G oc (iuii + iv n ). 

Inserting this into the Eq. (B21) yields the final form 
for the lattice susceptibility as stated in Eq. (15). 

As stated in the main text, the Matsubara-matrix 
formalism is not used in this work. Instead, the an- 
alytic continuation of all Matsubara frequencies to the 
real axis is the principal goal of the above treatment. 
This seems possible, as the local susceptibility function 
X loc («i'n|jwi, 10J2) along with its analytical continuations 

to the real axis x loc (^Vi ± z0 + |wi ± i0 + ,cj2 ± *0 + ) could 

be obtained from the effective impurity model. Equa- 
tion (B19), along with the irreducible scattering ma- 
trix (B22), can be continued to the real axis and the 
integral equation can be solved for the local two-particle 
Green functions £ 2 - A (v n ± i0 + |wi ± i0+ , co 2 ± i0 + ) . These 

are then used to determine the lattice susceptibility func- 
tion at the real axis via Eqs. (B13) and (B12). 

In practice, this procedure is rather involved and more 
importantly the extraction of the analytically continued 
two-particle Green function x loc (V„ ± «0 + |wi ± i0 + , ui-2 ± 

i0 + ) from the impurity model is far from trivial and 
strongly depends on the impurity solver. It is straight 
forward (but very cumbersome) within semi-analytical 
approaches like ENCA, where one has a direct handle on 
the analytic structure of the this function. But, for ex- 
ample, in QMC or NRG it is conceptual not clear how 
extract two-particle quantities at the real axis. 

In order to circumvent this involved procedure, we re- 
vert to an additional approximation. The summations 
over inner Matsubara frequencies are decoupled in a man- 
ner similar to the random phase approximation (RPA), 



^2 A{iv n I iu>3,i^2)Miu n | iu>i,iuj 3 ) (B27) 

iu>3 

=> Miv n I icJi,iu 2 ) y^^RjiVn |^i,^3)- 



TT = - IT&T) 



1 + n (t <g> r) 



(B24) 



With this, the inner frequency summations over iwi 
and iu>2 can be performed and the equations can be solved 
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cxplicity to yield the final form for the susceptibility dis- 
played in Eq. (16). 

The local and lattice particle-hole propagators enter- 
ing this equation can be calculated at the real axis us- 
ing standard techniques for the evaluation of Matsubara 
sums, 



£ loc (z/) = du f{u) p(uj)®G{u+v+i5) (B28) 



+ G(cj — u—iS) <g> p(uj) 



E(q,> 



oo BZ 



d D k 



(B29) 

p(k, uj) (g> G(k + q, uj + v + i8) 



(2ir) D 

+ G(k,uj — v~i8) ® p(k + q,uj) , 



with the spectral functions 



p(k, w) 
p{u) 



1 

d D k 



BZ 



(2tt 



(G(k,cj+iS)-G(k,w-iS)) (B30) 
D p{k,u). (B31) 



All momentum integrals are over the whole Brillouin- 
Zone (BZ) in D space dimensions. 
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